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ABSTRACT 



A method is described for determination of stresses 
near the circumferential junction of shell, head, and skirt 
of a pressure vessel subjected to mechanical and thermal 
loadings «, Considerable latitude is permitted for variation 
in geometrical, mechanical, and loading characteristics; 
however, it is presumed that there is complete symmetry 
about the axis of the vesselo A digital computer program 
is appended to facilitate solution of the problem.. 
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1. Introduction 

The chemical and petroleum industries have constructed 
many pressure vessels with vertical axes and cylindrical 
skirt supports. Some of these vessels contain fluids at 
temperatures as high as 1000°F with resultant high heat 
flows through the skirt support .Q0 V These heat flows, and 
the thermal stresses resulting, give rise to the "Hot 
Skirt Problem" which is here generalized to that of deter- 
mining the stress distribution near the junction of shell, 
head, and supporting skirt of a vertical pressure vessel, 
resulting from any combination of the following loadings? 

1. Changes in temperature proceeding axially along 
any of these elements (shell, head, skirt). 

2. Changes in temperature (assumed to vary linearly) 
through the thickness of an element. 

3 o Internal pressure in the shell and head but not 

in the skirt. 

* 

4. An additional axial loading on any or all of the 
elements. (Could result from the dead weight of 
one of the elements . ) 

5. A moment applied to the junction of the three 
elements by agencies external to the elements. 

6. A radially outward distributed loading applied to 
the junction by agencies external to the elements. 

^Numbers in brackets refer to bibliography, page 40. 
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Throughout it is assumed that the three elements are 
joined on a common circumference, that they and their load- 
ings are axially symmetric , that the material of the 
elements remains elastic, and that the structure always 
remains stable and all deformations are " small " „ It is 
also assumed that the walls of the elements are "thin",, 

The theory of thin axisymmetric cylinders (Section 2) 
shows that all except axial loadings applied a sufficient 
distance from the end of the cylinder (i.e», junction end), 
influence that end very slightly*, Axial loadings however 
are transmitted without change,, Also one may consider 
such shell configurations as cones, ellipsoids, etc«, as if 
they were cylinders for purposes of calculating their 
reactions close to the junction end due to influences at 
that end, providing such shells are reasonably "deep” o' 1 " 

Only a finite length of the shell need be considered since 
any nonaxial loading has negligible effect upon the junction 
end if it is remote enough „ This nonaxial loading may 
therefore be taken to be sero, which suggests a finite 
length with the remote end free,, 

Thus the present solution involves three finite 
cylinders, joined together at their common junction ends,. 

To be as general as possible, this solution allows the 
material properties of the shells to vary axially and even 

#See footnote, page 127, Ref„ 3 Part II„ 
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the radius and thickness to vary "slowly” , The variance 
of radius and thickness must be kept within reasonable 
limits to assure the desired accuracy in the solutions, 
since terms resulting from axial derivatives of these 
quantities are not included. 

Each cylinder is dealt with in turn, performing the 
following sequence of calculations: 

1. The moment, shear, slope (belling), and radial 
deflection are found and recorded for each of 
n (a preselected integer) points along the 
disconnected finite cylinder, resulting from the 
application of Loads 1 through 4 as previously 
listed, 

2. The same quantities are calculated and recorded, 
resulting from a unit moment (unit moment per unit 
circumference) applied to the junction end of the 
cylinder, 

3. The same quantities are calculated and recorded, 
resulting from a unit shear (unit radial load per 
unit circumference) applied to the junction end 
of the cylinder. 

Equations of statics may be written relating the end 
moments and shears to the externally applied loads numbers 
5 and 6, Using the junction end slopes and radial 
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deformations found in Calculations 1, 2, and 3, equations 
of slope and displacement continuity may be written 
relating again the end moments and shears <> These relations 
result in six linear equations with six unknowns (the 
three end moments and three end shears) which may be solved 
for the junction end moment and shear on each cylinder „ 

When the end conditions are known for an element , the 
deflection s moment, and shear along the length of the 
element may be calculated using the results of calculations 
1, 2, and 3, which have been retained . The stresses are 
then calculated using the deflection and memento 

Two methods are presented for handling internal 
pressure. The first considers pressure loading at the 
outset as an input into Calculation 1. This is permissible 
if the local geometry is essentially that of a cylinder. 

If, however, either the shell or head is not cylindrical 
in shape at the junction, all loads but pressure are 
included in Calculation 1. The effects of pressure are 
computed separately for shell and head using membrane 
theory,, The resultant end deflections due to pressure are 
included in the continuity equations and the membrane 
pressure stresses are added to the stresses calculated from 
all other loadings (including the bending stress caused 
by dissimiliar deflections due to pressure). 
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A finite difference iterative procedure is used to 
solve the fourth order differential eauation presented in 
Section 3 and is the "heart” of the solution 0 Essentially 
an assumption of the radial deformation of the cylinder is 
made and an integration is performed to find a new shape „ 

To assist convergence, an averaging process is used but for 
some shells convergence may still be difficult to attain,. 

The original goal was to develop a method which 
could be easily performed on a desk calculator in a 
reasonable length of time„ Due to the difficulty in 
attaining convergence this goal has only been partially met„ 
If one had a problem with cylinders of constant parameters 
so that the numerical solution need only be performed for 
Calculation 1, with an analytical solution available for 
Calculations 2 and 3, a hand solution might be feasible „ To 
solve more complicated problems a digital computer program 
has been written in Fortran language and is listed in 
Appendix B„ This program allows one to consider shells and 
loadings, the solutions for which would take a prohibitive 
amount of time by hand„ 

Three cylinders were assumed to be present in developing 
the method and the computer program,, However if one wishes 
to consider a vessel of but two abutting elements, one 
may choose a very small length for the unwanted cylinder. 
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thereby effectively eliminating its influence. 

Some previous work has been published on a hot skirt 
problem. Weil and Murphy £lj have developed a solution 
adapting the equations of the analytical solution of a beam 
on an elastic foundation to a cylinder. However these 
equations require a constant "beam" or constant cross 
section cylinder. This restricts the problem to a constant 
axial thermal gradient, and constant temperature difference 
across the wall thickness. These restrictions are severe 
and limit greatly the usefulness of their solution. 
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2. Theory 

Consider a small element cut from a cylinder with 
dimensions dx, ad , and t in the x, y, and z directions . 
The cylinder is of radius a, and has a wall thickness t, 
Figure 1 illustrates the element . 




Figure 1* Element Cut from Cylinder 
The following symbols will be used in this section: 
a - radius of cylinder 

D s Et^/12 (l-/^ ) = longitudinal stiffness of a unit width of 
element, (also called plate flexural rigidity „ ) 
d = infinitesimal central angle 0 
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E - Young T s Modulus c 

and e„,_ s centerline extensional strain in x and z 
mx. mz 

directions respectively due to all but thermal 
causes . 

e bx and e bz - bending strain on inner wall surface in x 
and z directions respectively due to all 
but thermal causes; and positive as indicated 
in equations (5) and (6). 

e x and e z - extensional strain in x and z directions 
respectively., 

M x = moment (as illustrated) per unit circumference . 

M z - moment (as illustrated) per unit axial length. 

P = ap/2 + P' ,s s axial force per unit circumference 0 

P v = axial force per unit circumference from all causes but 
pressure. 

p - internal pressure. 

s - radial deflection. (Positive outward.) 

t - wall thickness. 

T - temperature at centerline of wall. 

AT z inner wall temperature minus outer wall temperature. 

V = radial force per unit circumference. (Positive as shown.) 

x, y, and z = cartesian coordinates as shown in Figure 1. 

o< z coefficient of linear thermal expansion, 
s curvature in x, y plane. 
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p x z radius of curvature in x, y plane, 

C T x and 0~ z s axial and circumferential tensile stresses, 
respectively, 
jjl = Poisson^ ratio. 



For axial equilibrium: 
't/2 



■X 



(1) P 

For equilibrium in the y direction: 



/ 

t/2 



pad£dx +■ ad0V = ad0(V+-d?) + d0dx 



/ 



t/2 

0“ dy 
t/2 2 



This reduces to: 




Taking moments about axis o~o: 



ad£(M x + Vdx) * (pad0dx)“ s (N^ + dM x )ad0 £ 

(d#dX /t/2 <I ' zdy) f 

which reduces to: 



(3) 




dx 



Also from elementary statics : 



(4a, b) M x 




CT x ydy and M z s 




0~ z Y&y° 
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Assuming the longitudinal strain is linear through 
the wall: 

( 5 ) s ©rnx ” “^bx “ ^ T + ©CT 

t t 

Assuming circumferential strain is linear through 
the wall: 

(6) e z = e mz - “^bz “ ^T 4 . oct 

T> t 

These assumptions are reasonable if t<S£a. Otherwise 
the distribution may be nonlinear . Thus ? for example s the 
stress distribution due to internal pressure in a thick 
cylinder is hyperbolic (Lame / Solution). 

Figure 2 shows the cross section of a piece cut from 
the cylinder wall of width ad£ and length dx. The strains 
are greatly magnified for illustrative purposes. From 
similar sectors: 

X x dx s dx /px s' ( e bx + otAT/2)dx/(t/2) . 
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Figure 2 Cross Section of Cylinder Wall 

Presuming that ii^^l, Therefore: 

dx * dx 2 

(7) 3% _ ( 2e +*4T)/t. 

dx*~ bx 

(S) e 2 = s/(a+y)^s/a. (Here again we assume t^a„) 

Comparing (6) and (S), and equating coefficients of like 
powers of y: 

(9) e bz - -~AT/2 

(10) e m2 = s/a - •<T. 

The fundamental stress-strain relations are: 
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(11) £e/(1-/c )] [j e mx' < ’^ e mz^ ” ^f^ e bx^ /* e bz iJ 

(12) C^z” [E/(l-/t)J Q e m z +/<. e rDx^ " e bz^/^ e bx 0 

Therefore: 



(13) 


P = 


A/2 

/t/2^ dy 


= [et/C 


(14) 


dV 

dx 


“ P 


• *J 


’t/2 

.t/f 2dy 


(15) 




/ 


'%/ 2 




Mx 


■-1 


, cr ydy = L Et 

■t/2 X 


(16) 


M 

z 


-i 


't/2 r 

.t /2 ^ ydy = iEt 




Combining (7), 


(9), and 


(17) 


d^s 

dx^ 


~ D 


(l+^Kd T 
* t 



mz v~" --mx 



‘i--” 2 !] H e bx*/< e bz , J 



Eliminating e M from (13) and (14): 



(18) g »p 



Ete 

/<r mz 



a a 

Substituting (10): 

dV „ /*P . Etoc T 

(19) ax = p - ~ * ~T~ 

From (15) and (9): 

(20) e bx s 6( 1 )M * + >**• A 



Ets 

„2 



Et 
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and from (13) and (10) 



(21) e s 11? ^1~L- - /4£s 

mx Et a 



Summarizing: 

, x dV _ „ A Et<x.T 

(22a) = P + —t— 



Ml 

a 



Ets 

_2 



(22b) s 
dx 



(22c) 



d 2 s 

dx 2 



% + (1-^)X4T 

D t 



If it is assumed that T-P- AT“0 S and that a, E s/ # , and 
t are constants, equations (22) reduce to: 



(23) 

dx 1 * 



>Ets 



is : 



The solution to (23) assuming a semi-infinite cylinder 



fix 



Et 

4a 2 D 



(24) s - —j— V^cos^sx ^^R^jtcos^x - sin^x) 

2^ 

where and are the end loadings, and ^ s 

d 2 s 

M„ then equals D-— e 
x dx 2 

After Mjj. and s have been determined either by numerical 
solution of (22) or by using (24), the stresses in the 
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and from (13) and (10): 



(21) e s . (l-// 2 ) ,P _. - 

mx Et a 



Summarizing: 

, dV _ „ A EtocT 

(22a) gj = p ♦ — r- 



M 

a 



Ets 

„2 



(22b) s 
dx 



(22c) 



d 2 s 

dx 2 



% + (l-^)^^T 



D 



If it is assumed that T^Ps and that a, E s/ #, and 

t are constants, equations (22) reduce to: 



(23) 

dx 4 



-Ets 



is : 



The solution to (23) assuming a semi-infinite cylinder , 



1 fi x 



(24) s - — j — V^cos/Sx f^5M xl (cos^gx - sin^x) 

where Vl and are the end loadings, and ^ s 



Et__ 

4a 2 D 



M x then equals D : 



d 2 s 
dx 2 



After and s have been determined either by numerical 
solution of (22) or by using (24), the stresses in the 
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cylinder may be determined with equations (9) , (10) , (20) , 
(21), (11), and (12). The results are: 



CT 



x inner = 



P 6M x 
t .2 



(25) 



0 “ 



x outer = T ~ 



6M 



CT" . Es ,/a P 

z inner = — E*T + / ^ 



( 26 ) 



Es 

a 



outer = — ~ ~ E*T + '’ir 



M? 

t 



2 t 2 



Ep<4T 
+ 2 




o 



A single equivalent stress is frequently used as a 

basis of assessing the severity of loading on a structure 

in which the state of stress is not a simple uniaxial 

stress. The equivalent stress employed here is- that of 

the distortion energy theory (equivalently, of the octa~ 

hedral shearing stress theory) . If the principal stresses 

are (T~ , 0~ and 2 ero, as they are at the outside of the 
x 2 

cylinders, the equivalent stress is: 




On the interior surface, if a pressure p is acting, the 

third principal stress is ~p rather than zero. However, 

the same formula is used since CT and 0~ are at least one 

x z 

or two orders of magnitude greater than p. Thus: 
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(27) 

and : 
( 23 ) 




outer 



inner 



(7** - (£ TO") +< T c 

x outer x z outer z 






- (C TCT). +(T fc 

x inner x z inner z 
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3. 



Method 



Figure 3 shows the basic three-cylinder problem that 
is to be solved, with the numbering convention used, and 
the sign conventions for the externally applied moment, 
shear force, radial deflection, and slope when the cylinders 
are considered jointly. 



‘t Positive radial deflection for all cylinders. 




Positive slope for all cylinders. 




Positive shear for all 
cylinders. 

M Positive moment for all 
cylinders. 

Cylinder 3 



+x _ ax is 



0 



_ a xis _ +x 



Figure 3. Positions of Three Cylinders 
and Sign Conventions 
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The basic procedure to be followed in the solution is 
to solve equations (22) numerically or (23) analytically 
to find the reaction of each cylinder to (1) the given 
loadings, (2) a unit junction end moment#, and (3) a unit 
junction end shear# when this cylinder is "free" (i°e°, 
disconnected from the junction) « The three cylinders are 
then conceived to be rejoined, and by using the above in- 
formation along with the equations of statics and continuity, 
the reactions (i°e°, the actual end loadings) of the three 
cylinders are determined „ 

The analytic solution of equation (23) is given in 
equation (24)° However in the numerical solution of (22) 
we choose a finite length for the cylinder,. The only load 
applied to the remote end is a moment equal to that caused 
by the radial thermal gradient. From equation ( 24 ) 

the radial deflection at x = 5/p is about 0°7$ of the 
radial deflection at x = 0 for any given loadings applied 
at x = 0° Thus, if L, the length of the finite cylinder, 
is set equal to 5/p there is only a small error in the 
radial deflection at the junction end due to neglecting 
the loading applied at x = 5/p by the portion of the cylinder 

""^Throughout this thesis "shear" means the radially 
outward directed load per unit circumferential length and 
"moment" means moment per unit circumferential (or axial) 
length* 
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beyond this point. However an error results if there are 
nonlinearities close to the remote end. See ’’Discussion 
of Results,” Appendix C. 




^ 0 axis 




Figure 4, Finite Cylinder Divided into n=l 

Increments Showing Loads at Junction 
and Remote Ends. 

The following iterative procedure was developed by 
John Eo Brock, Professor of Mechanical Engineering, U. S. 
Naval Postgraduate School, Monterey, Calif., and modified 
by the writer to assist convergence. 

Figure 4 shows a cylinder divided into n~l equal 
increments with stations 1 through n. (With a modification 
of the integration method used, the requirement of equal 
increments could be withdrawn.) At each station, values 



IS 



of all the parameters of the cylinder and all loadings are 
assumed known . 

Equations (22) may be restated as follows: 

+ EtyxT //P Ets 



(29) 



dV 

dx 



a 



a 



v =, 


f dV 
• cE dx 


M = 

X 


jVdx 


d 2 s 

s? 


M 

X 

= t r + ■ 


ds . 
dx 


f d2 fdx 
y dx 



Li '/ix 1 * 

Let M* = M + 

x t 

(30) M* = /Vdx + D(1+/0*4T 

* t 

(31) d£s = M* 

,2 D 
dx 

Let s = A + Bx + f(x) where A and B are constants and 

f(x) is normalized in such a way that f(0) = f(L) = 0 o 

Also let Q = o Then from (29) , (30), and (31): 
dx 

Q = p + _ / ^P_Etj7 + Bx + f( x j] 

a. 3 - ^ 



(32) = Q - A% - . Etfixl 

a or a^ 
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where Q = p + 

a a 

✓ 


fX-^ 


V = (Qdx = V +J Qdx - Aj 


1 dx “ 




> a 


Let ^ = k 

a 4 

/* r x 




V = V - A/ kdx - B / kxdx 


- / kf (x) 


_ /0 fX _ 

where V = Qdx 0 


Jo 


Next 





M* = M + j 0 Vdx - A^fjkdxdx - By^fkxdxdx 

D(1+mMT 



dxdx 



so that 



(34) M* = - Aj^^kdxdx - Bj^^kxdxdx J 5 kf(x)dxdx 

where M* = M + ( Vdx + 

xl Jo t 

An initial assumption is made for f (x) ; conveniently, 
one may take f(x)=0. Therefore the only unknowns in 
equations (33) and (34) are A and Bo 

It is also assumed that at x = L (the remote end of the 

cylinder) , M* = V n = 0» This is a reasonable assumption 

because as x approaches infinity, with a constant radial 

thermal gradient, the axial curvature of the cylinder wall 

approaches zero» L has been chosen long enough so that the 

conditions at x = L approach those at infinity which allows 

M* = V = 0» The two equations and two unknowns are: 
n n 
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(35) 0 



= V n ” A X kdx “ B L kxdx " X 

-*■ /X 



L 



o 

(36) 0 = Mh - A X X kdxdx - B J 0 kxdxdx = j t ^kf(x)dx 
Solving for A and B: 

[V n - y^kf (x)dxj [i /kxdxdxj - \J 0 kxdxj Jkf ( x) dxdx ' 



kf (x) dx 



(37) A = 



(3d) B = 



[Jj kdxj //C Jo kxdxdx/ “ jjkxdx] [j » X kdxdxj 

[f 0 kdx] fe-i i''kf(x)dxd^--[v n ^kf(x)dx/[X^kdxdx 



KH [ JjJo kxdxdx j kxdx j fXt kdxdx 7 

Substituting these values into (34) gives M*. 



(3D 



d 2 s _ M* 

. 2 * D 
dx 

Ij 0 t O 



(39) g(x) 



. A 2 s 



’£ 



dx = 



dx 

* 



_/V 

'A ~b 



dx 



(40) f-^x) =-j Q g(x) dx 

It will be noted that g(x) ^ ^ and f^(x) / s because 

the initial conditions have not been included . We normalize 
f^(x) to obtain 

( 41 ) fjx) = f _ ( x) - f l< L ) x 



and it may be seen that ££(0) = f 2 (L) = 0. 



If f£( x ) is used for f(x) directly in equation (32) to 
perform the next iteration, the process usually diverges 
rapidly. To force convergence, the following averaging 
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procedure is used 



lo Let numbers m and C . . chosen as described later, 

mm* 

be specified., 

2. Let d be the value of x for which f.^Cx) is maximal . 
3o Determine a number C according to the rule: 



mxn s 



If f(d) = 0, C = C • : 

' ’ mm’ 

If |f(d)/f 2 (d)| > 1, C = C r 

If 0<|f ( d ) /f ^ C d ) | < 1 s C equals the larger 

of C min or I r(d)/r 2 ( d ) | ». 

Obtain a new function f^(x) by the formula: 



(42) f^ (x) = Cf 2 (x) + (l-C)f(x)o 

f^(x) is the function to be substituted into the equation 
(32) for f(x) in the next iteration., The process has 
converged when f 2 ( x) (x) . 



The fraction 



f(d) l m 
?iTdT| 



is used solely to reduce the 



number of iterations required for convergence . It may be 
seen that when convergence has almost been attained, the 
fraction C will become progressively larger resulting in 
a higher percentage of f 2 (x) being used in the next 
iteration. 

The number m must be kept to a minimum and C . kept 

mm r 

to a rraximum to reduce the number of iterations to a 
minimum. It was found by trial and error that a value of 
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6 or 9 for m and <,06 to „ 0 S for C min generally insured 
convergence o One problem required an m of 15 for reasons 
unknown to the writer „ Generally if divergence is exper- 
ienced, increasing m and decreasing C m ^ n seems to eliminate 
the divergence . 

When convergence has been attained, M , s, and 4§. are 

x 

calculated as follows: 

(43) M x = M*- DU- 1 * NIT 

(44) gf = B + g(x)- 



(45) s = A +Bx + f 2 ( oc) 

Using the sign conventions established in Figure 3 and 
with the understanding that ’’loadings" stand for the pressure, 
thermal loadings, and axial loading in each cylinder, let; 
< ^ol , ®o 2 , ^o 3 = en d deflections of the disconnected 

Cylinders 1,2, and 3 respectively due to their 
loadings . 

^1 ,^2 = junction end deflections of Cylinders 1,2, and 

3 respectively caused by joining the cylinders . 
at the junction o The total deflection of 
Cylinder 1 would then be « 

4 > ol , ^>3 = j unct ^ on en< ^ slopes (subscripts have the same 

$1 , ^ > ^3 meanings as those for £ ) « 

M 0 i,M 0 2 >M 0 3 = junction end moments of Cylinders 1,2, and 3 
respectively* 
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V 0 ijV 0 2 >V 0 3 = junction end shear loadings of Cylinders 1,2, 
and 3 respectively 0 

V 0O and = externally applied shear and moment at 
junction o 

= slope at junction end of Cylinder 1 produced 
by application of unit shear to junction ends 
A v £ = deflection at junction end of Cylinder 1 
produced by application of unit shear to 
junction end* 

= slope at junction end of Cylinder 1 produced 
by application of unit moment to junction ends 
A m g = deflection at junction end of Cylinder 1 
produced by application of unit moment to 
junction ends 

B and C denote similar coefficients for Cylinders 2 and 3 
respectively, with the same subscripts and definitions as 
for A. 

Mote that <P oV <t> oZ , k vi , A mr B^, B m( ,, C mlM and C m f 
are the negatives of the end conditions calculated by the 
above iterative process, or analytical process, due to the 
difference in sign conventions between Figures 3 and 4° 
Figure 5 shows Cylinders 1,2, and 3, situated as in Figure 
4 but with sign conventions as in Figure 3, and as used in 
these equations o 
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y 




A 



Positive Positive 

0 Slop e \ Defle ction a xis x 



y 




Figure 5. View of Cylinders 1,2 , and 3> Showing 

Sign Conventions Used in Equations (46), 
(47), and (4S)» 
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From statics, compatibility, and other considerations 



(46) 



V ol + V o2 + V o3 



00 



+ M 0 +M=M 
ol o2 o3 oo 

+ Ki - K + Kz 

^1 + 4l = S + ^o3 

*1 + *ol = ^2 + *02 
*1 + *ol ” *3 + *o3 
^ol A W> + ^ol A m3 “^L 
V ol A vS + M ol A m£ “^1 
V o2 B v«. + M o2 B m* ~ S Z 
V o2 B vS + H oZ B mi 
T o3 C v< + M o3 C m^ = *3 



V 0 C . + M .C f 

o3 v£ o3 m£ 






Solution of equations (46) yield: 



M 



_ p 2^1 ~ p 1^2 



o2 



N 1 P 2 



P 1 N 2 



M 



o3 



= N lQ2 “ N 2^1 



(47) 



N P 
12 



N P 
2 1 



M = M 
ol oo 



M 



M 



o2 o3 

v o3 - U/ipUi - JiM o2 - k A3 ) 

V o2 - (l/DiXHi-E^ - F lMo2 - G lMo3 ) 



V . « V 

Ol OO 



o2 



V 



o3 



where 
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D 1 


= 


A v4> + 










D 2 


= 


A v4> 










d 3 


= 


A vS + 


B vS 








d 4 


= 


A vS 










E 1 


= 


A v<P 










E 2 


= 


A v4> + 


c v <i> 








E 3 


= 


A v$ 










E 4 


= 


A vS + 


C v< 








F 1 


= 


A m <t> + 


B m d> 








F 2 


= 


A m<^ 










F 3 


3= 


A m S + 


B m 6 








F 4 


= 


A m 6 










G 1 


= 


A m<i> 










°2 


= 


A m <*> + 


^m<t> 








°3 


= 












G 4 


= 


A mg + 


G , 
mg 








H 1 


= 


*01- 


*02 + 


V A „ 
00 V P 


+ 


M A 
oo m<* 


H 2 


= 


(pd 


*03 + 


V A . 

00 V® 


+ 


M A j 
oo mi) 


H 3 


= 


8 ol- 


S 02 + 


V A 
oo v6 


+ 


M A 
oo m S 


H 4 


ss 


8 ol 


S , + 
o3 


V A 
oo v£ 


+ 


M A 
oo mf 


h 


= 


D 2 E 1 ' 


- D 1 E 2 








h 


= 


D 3 E 2 ' 


' D 2 E 3 








*3 


= 


D 4 E 3 ' 


‘ D 3 E 4 
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(48) 



J 1 


— 


Vl 


— 


D F 
12 


J 2 


= 


D 3 F 2 


- 


D F 
2 3 


J 3 


= 


D 

4 3 


- 


D 3 F 4 


K 1 


= 


D 2 G 1 


- 


D 1 G 2 


K 2 




D 3°2 


- 


D 2 G 3 


K 3 


= 


D 4°3 


- 


D 3 G 4 


L 1 


= 


D 2 H 1 


- 


D 1 H 2 


L 2 


= 


D 3 H 2 


- 


D 2 H 3 


L 3 


= 


D 4 H 3 


- 


D 3 H 4 


H 1 


= 


Vl 


« 


I 1 J 2 


N 2 


= 


V 2 


- 


I 2 J 3 


P 1 


= 


Vi 


- 


X 1 K 2 


P 2 


= 


Vz 


- 


I 2 K 3 


«1 


= 


Vl 


- 


I 1 L 2 


«2 


= 


Vz 


- 


X 2 L 3 


Let 


■ m 1x- m 2x and M : 



Cylinders 1,2, and 3 respectively. 



M L1* M L2 * 



and M. 



L3 



represent the axial moments in Cylinders 1,2, and 3 caused 
by the thermal, pressure and axial loadings . M M V 2 8 
and represent the axial moments in Cylinders 1,2, and 
3 respectively caused by application of a unit shear load 



to their junction ends» M^, M m2 , anc * M m3 re P resent the 

axial moments in Cylinders 1,2, and 3 respectively caused 
by application of a unit moment to their junction ends* 
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The subscripts for V and s have comparable meanings. 
Then: 

M lx = M L1 + V ol M vl + 

K 2x = M L2 + V o2\2 + M o2 K m2 



(49) 



M 



3x 



M- 



l L3 + V o3 M v3 






■3 M m3 

V lx - V L1 + Wvl + M ol V ml 
V 2x * V L2 + V o2 V v2 + M o2 V m2 
v 3x = V L3 + V o3 V v3 * M o3 V m3 
s lx “ S L1 + V ol s vl + M ol s ml 
s 2x = S L2 + V o2 s v2 + M o2 s m2 



s 3x " S L3 + V o3 s v3 



M^-,s 



o3 m3 



Substituting the results of (49) into (25)? (26\ (27) 9 and 
(2S) gives the desired stresses along the cylinders » 

In the above development it was assumed that the head 
and shell (Cylinders 2 and 3) were in fact cylinders and 
that their pressure stresses could be computed accordingly., 
Call this method Option 1* If one or both are not cylinders 
but some other shell of revolution (ellipsoid, cone, etc.) , 
an error will result in the stress computations by consid= 
ering them cylinders . (See problems 10 and 20, Appendix Co) 
One way to reduce this error is to use the above methods 
to obtain the reactions resulting from all loadings but 
pressure, and then superpose the separately calculated 
reactions due to pressure. Call this method Option 2. 
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Figure 6. Shell Head, and Skirt Showing R1 and R2. 

Consider the shell and head shown in Figure 6. In 
this particular case the head is Cylinder 2 and the vessel 
shell is Cylinder 3 but the following development will 
allow a noncylindrical element for either the shell, head, 
or both. The vessel shell is a cylinder, the head is some 
other shell of revolution. R is the radius of the vessel 
at the junction, R^ is the radius of curvature of the 
meridian of the head at the junction, and R^ is the radius 
of curvature of the line in the head perpendicular to the 
meridian at the junction. From Figure 6, -- Cos 6*, where 

0 is formed by the axis and the tangent to the meridian at 
the junction. 
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From membrane theory: [ 3 ] 



( 50 ) &2 _ £ 

R 1 R 2 ‘ " 

where ^ and ^ are the meridional and circumferential (or 
hoop) stresses respectively, p is the inner gage pressure, 
and t is the wall thickness 0 At the junction the axial 
component of ^2 must equal to satisfy staticso There- 
fore : 

pRR 2 pR 2 



(51) ^ = 



EL 



2tCos0 2tR 
From ( 50 ) : 

(Ti 



2t 



(51a) G T 2 - R 2 [§ -^= R, [§ - If] 

Let be the radial deflection due to pre 



ssure , 



<52) ^-[lOU 



E Ai£] 

2tR 2tJ 



There is also an inward radial force at the junction, 
due to the component ofG^ perpendicular to the axiso 



(53) V£ q =-0^tSin©= -p^ 2 “ R ] 



The suggested procedure for computing the effects of 
pressure in Option 2 is first to calculate the inputs to 
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equations (48) setting internal pressure equal to zero but 

considering all other loadings and external loadingSo £ ® 

o2 

and <C® are then calculated and added to £ „ and S , o V® 

^o3 o2 o3 oo 

is calculated as the sum of the contributions of both 

Cylinders 2 and 3 and is added to V 0 Using equations 

U V 

(47), (49) 9 ( 25 ), and (26), the longitudinal and circum- 
ferential stresses are computed,* However, it will be noted 
that these stresses contain only those pressure stresses 
caused by the dissimilar deflections and consequent bending 
of Cylinders 2 and 3° The direct membrane pressure stresses 
must still be added as follows; 





^x 


outer = 




outer 


(from 


(251) 


+ r i 


(54) 


<Tx 


inner = 


<r* 


inner 


(from 


(25)) 


+ 0"i 




outer ~ 




outer 


(from 


(26)) 


+ <r 2 






inner " 


% 


inner 


(from 


(26)) 


+ 0 ~ 2 



These values may then be substituted into (27) and (28) to 
obtain the equivalent stresses*. 

Note that with Option 1 the effects of internal press- 
ure are calculated in the numerical solution of the diffei*- 
ential equations and therefore any variations in the 
parameters of either cylinder will be taken into account „ 
These variations will cause moments and shears in the 
cylinder itself before consideration is given to its junction 
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with the others o With Option 2, the only moments and 
shears due to pressure that are calculated are those caused 
by the dissimilar deflections at the junction end of each 
cylinder., Either option treats all three shells as cylinders 
when calculating the effects of other than pressure loadings.. 
For this reason either option will err in treating axial 
force in a shell with other than zero meridianal curvature 
since, among other things, the cylindrical equations do not 
consider the radial component of CT Xo 
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Digital Computer Program 

At first the goal of this thesis was to develop a 
method in which a complicated hot skirt problem could be 
solved on a desk calculator in no more than a day 0 This 
goal was not met due to the number of iterations necessary 
for convergence o 

To make the method more useful 8 a digital computer 
program which closely parallels the method outlined in the 
previous section was written in the GDC 1604 Fortran 60 
language o There are some small deviations in sequence of 
operations mainly due to the fact that the program progress- 
ively grew as it became more apparent that a manual solution 
was becoming impractical for any problem in which the 
parameters of the cylinders varied axially <> A standard 
library program is used to solve equations (46) » 

The program also provides for longitudinal slots in a 
portion of a cylinder (usually the skirt) „ Cheng and Weil 
[2] have indicated that slotting the skirt at the junction 
substantially reduces the stresses in the cylinders., 

The slots are considered narrow enough so that the 
loss of metal in them may be neglected 0 The slotted section 
is considered to be a series of separate cantilever plates 
with no interaction between adjacent plates and* conse- 
quently, no (T along their slotted lengths, where T ms is 
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the stress at the midplane in the z direction There is 
of course some at the junction end of a slotted skirt 

due to the strain of the junction itself and therefore 
this is not an exact representation of the actual situation u 
However provision is made in the program to start and end 
the slot at any distance from the junction one might choose, 
and a way might be found to simulate the junction conditions 
by starting the slot at some predictable distance from the 
junction., This slotted shell capability was provided 
mainly so that others may have a useful tool for system- 
atically investigating the effects of such slots.. 

The following alterations are made to the equations 
in the program at those stations where the cylinder is 
slotted: 

I. Q of equation (32) is set equal to zero since all 
terms except p are functions of interactions between 
adjacent longitudinal segments of the cylinder » If 
the cylinder were slotted it could carry no internal 
pressure 0 

2o For the same reasons all the terms except the 
6Mx/t/ ^2 term are set equal to zero in equations 

(26) o 

The program uses a midordinate integration scheme 
which, as programmed, wastes about half as much memory space 
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in the computor as it employs . It was originally used 
because it is systematic and very convenient for hand 
calculations o If limitations should thereby be imposed 
upon the utility of the program, more compact integrating 
methods could be used. 

The program was written so that only the minimum 
number of input cards necessary to describe the problem 
is required . For example if only the first and last 
stations of a cylinder are described, the program uses the 
cylinder parameters of the first station for the remaining 
intermediate stations and makes a linear interpolation of 
the thermal and axial load inputs. The internal pressure 
is assumed constant. See Appendix A for a detailed 
explanation. 

Three output options are available to the user: 

1. Only the maximum equivalent stress (based on the 
distortion energy theory) and its location is given 
for each cylinder. 

2. The inner and outer longitudinal, circumferential, 
and equivalent stresses are given for each station 
of each cylinder. 

3. The inner and outer longitudinal and circumferential 
stresses and the moment and shear are given for 
each station of each cylinder. 
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With each option all of the input data is repeated, and a 
short explanation of the output is printed <> 

The program takes about one minute to compile and 
each problem requires from one half minute to one minute » 
It took about seven minutes to compile and solve the 
eleven problems in Appendix Go The program uses only the 
normal Fortran functions and language as described by 
McCracken Cy and therefore should be easily adaptable to 
any large digital computer having a Fortran system „ 
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5 



Conclusion 



Even though the method was developed specifically to 
solve the "Hot Skirt" problem, the reader is reminded that 
with suitable manipulations of the input, many other 
problems may be solved., 

The only restrictions are the following 
lo Elastic behavior., 

2o Axial symmetry o 

3° Maximum of three shells arranged as in Figure 3® 

4o Constant internal pressure in Cylinders 2 and 3 only., 
5 o "Thin” wall* 

6o Linear temperature variation in radial direction,, 

7° Configurations of shells such that their reactions 
to an end load may be approximated by regarding 
them as cylinders of the same radius, wail thickness, 
and material o (Thus, treating a very shallow 
elliptical head as such probably would result in 
appreciable error «) 

Some useful extensions of the development would be: 
lo A complete revision of the theory so as to include 
a radius of curvature in the x.% plane 0 This would 
eliminate the restrictions imposed by the use of 
the cylindrical equations 0 
2o An alteration to allow inelastic behavior 0 This 



3 $ 



would probably take the form of a series of 
repetive stress calculations, increasing the 
loading until yielding begins, and then a different 
set of calculations until the full loading is 
applied® 

3® A modification to accommodate ’’thick" shells such 
as are found in reactor pressure vessels® 

4 . A modification to allow radially variable material 
parameters such as would be found in a reactor 
pressure vessel with internal cladding® 

5® Modification to include the moments due to axial 
forces acting through offsets in radii or through 
radial deflections® 
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APPENDIX A 



Instructions for use of Program Skirt 
These pages are intended to serve both as Appendix A 
of the thesis 

"Thermal Stress Analysis of Pressure Vessels 
with Cylindrical Skirt Supports" 
by James Po McManus, May 1964, 
and as an independent document describing the employment 
of Program Skirt to solve certain stress=analytical 
problems for pressure vessels 0 The theory upon which 
Program Skirt is based is developed in detail in the above 
mentioned the sis 0 Those interested in studying the larger 
document should make inquiry from Head Librarian: U 0 So 

Naval Postgraduate School, Monterey, California „ 

Given three axisymmetrie shells joined together and 
numbered as shown in Figure (Al)„ These shells may be 
subjected to the following axisymmetrie loads: 

lo Changes in temperature proceeding axially along any 
or all of the shells „ 

2 0 Changes in temperature (assumed to vary linearily) 
through the thickness of a shell „ 

3» Internal pressure in Shells 2 and 3 but not in Shell 

lo 

4° An additional axial loading on any or all of the 



A! 



shells o (Could be the dead weight of one of the 
shells and its contentSc) 

5. A moment applied to the junction of the three shells 
by agencies external to the shells <> 

6. A radially outward distributed loading applied to 
the junction by agencies external to the shells „ 



't- Positive radial deflection for all cylinders. 



Positive slope for all cylinders. 



Cylinder 1 




Positive shear for all 

UU „ 0 c, 

cylinders 



M Positive moment for 
all cylinders c 



Cylinder 2 



Cylinder 3 



±x axis Q-l 0 _ axis ±2t 

Figure Al„ Positions of Three Cylinders and 
Sign Conventions 

Program Skirt solves the problem using simple numerical 
techniques. The user must choose a finite length for each 
shell such that loads applied at the remote end of the shell 
have no appreciable effect on the junction end, except that 
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axial loads are fully transmitted 0 Suggested values for 
the lengths are given in the remarks on Card 2o Each 
finite shell is then divided into n=-l equal increments with 
the dividing stations numbered from 1 at the junction to n 
at the remote endo Data for each station such as shell 
parameters and loadings is either read into the computer, 
or created by the computer from data given about adjacent 
stationso 

Two options are available for computing the effects of 
internal pressure „ Either option may be used if the shells 
are uniform and are cylindrical or have small meridional 
curvature near the junction,, Option 1 is to be preferred 
if the shells are cylindrical or have small meridional 
curvature, but there are variations in thickness, material 
properties, etc. Option 2 is to be preferred if there are 
no variations of thickness or material properties but the 
meridional curvature at the junction is not small „ Neither 
option is capable of obtaining accurate results if the 
shell is not relatively deep„ 

The following assumptions have been made in the program 
1„ The three shells act as elastic cylinders „ The 
shells may have "slowly variable" radii, thickness, 
and material properties, but are still treated as 
cylinders from station to station „ 
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2« Temperature varies linearly through the wall thick- 
n 6 S3 o 

3« Shear and curvature are negligible at the remote 
endo If there are nonlinearities near the remote 
end such as variation in the shell parameters, non- 
linear axial thermal gradient, or a changing radial 
gradient or axial load, the values of stress 
obtained near the remote end will be inaccurate o 
Description of Input Cards 

The input consists of all the parameters necessary to 
describe the problem plus some control parameters to direct 
the method of solution and the format of the output,, In 
the descriptions below the items given from left to right 
are: Column Numbers, Format, Parameter, Remarks,, When a 

parameter has an ft I ,? Format Field, its last digit must be 
in the last column of its fieldo 

Card Type I (One card only) 

1-10 110 Problem number May be 1 through 999 „ If 

it is zero or blank, the 
program stops o 

11-20 110 Output option 1, 2, or 3 

number 1 Only maximum equivalent 

energy stresses 5 *' and loca= 



^Equivalent on basis of distortion energy theory „ 



21-30 


FlOoO 


Externally 
applied moment 
to junction 


31-40 


FlOoO 


Externally 
applied shear 
to junction 


41-50 


FlOoO 


Pressure 


51-60 


110 


Code for hand- 
ling pressure 



tion given for each 

cylinder o 

2 Longitudinal, circumfer- 
ential and equivalent 
stresses given for each 
station of each cylinder 

3 Longitudinal and circum- 
ferential stresses, and 
moment and shear given 
for each station of each 
cylinder o 

Positive if causing com- 
pression on outer surface 
of Cylinders 1 and 2* (Mo- 
ment per unit circumference » ) 
Positive if directed radially 
outward o (Shear load per 
unit circumference o ) 

Gage pressure in Cylinders 
2 and 3o 

If 1, Option 1 is use do (See 
above o) If 2, Option 2 is 
usedo 
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61-70 F10.0 Minimum 

fraction 



Card Type 2 



1- 5 


15 


Cylinder number 


6-10 


15 


Number of 






increments 


11-15 


15 


Interval of 



printout for 
input data 

16-20 15 Exponent for 

forcing 
convergence 



The minimum value of a frac- 
tion C . used in the 
mm 

averaging routine for assis- 
ting convergence o Recommend 
from oQ 6 to „08. If left 
blank, .06 substituted. 

Values less than o 06 might 
be necessary if divergence 
is experienced. 

(Three cards required) 

1, 2, or 3. Must be in 
sequence . 

Maximum of 100 . 

Example : If 2, every other 

station reported; if 5, every 
fifth station reported, etc. 
Suggest 8 to 12 o If diver- 
gence is experienced, use a 
higher value. The higher 
the exponent , the greater 
the computer time re quire do 
11 substituted if left blank 



A6 



21=30 FXOoO 
31=40 FlOoO 
41=50 FlQoO 

51=60 FlOoO 
61=70 FlOoO 



Minimum x 

Maximum x 

Convergence 

criterion 



Value of x at 
beginning of 
slots 

Value of x at 
end of slots 



See belowo# x at junction 
end* 

See belowo* x at remote 
end<, 

Suggested values from O 01 to 
o05« This is an approximate 
indication of the relative 
error in the results 0 
Smaller values require 
longer computer time 0 
For use with longitudinally 
slotted sylinderso 

For use with longitudinally 
slotted cylinders o 



Card Type 3 (Two or more cards for 

each cylinder) 

1= 5 15 Station number 1 through n, where n is the 

number of increments plus one » 



*The program allows x, the axial dimension in inches, 
to be other than zero at the junction if the user so desires « 
It is recommended that the length of the cylinder be approx- 
imately six times the square root of the radius times the 
wall thickness, or 6Vat 0 
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6-20 


F15-0 


Young * s Modulus 
at this station 


21-35 


F15.0 


Mean radius of 

cylinder 

at this station 


36-50 


F15-0 


Poisson’s ratio at 
this station 


51-65 


F15-0 


Wall thickness at 
this station 


66-S0 


F15o0 


Thermal coefficient 
of expansion at 
this station 



Note: There must be Cards 3 for the first and last stations- 
The first station is number 1 and the last is n, If 
there is no card for an intermediate station, the 
values for the previous station are substituted- 
Therefore to change a value, simply feed in a card 
for the appropriate station- The cards submitted 
must be in numerical sequence by station number- The 
minimum number of Cards 3 for a cylinder is two, the 
maximum number is n- 

Card Type 4 (Two or more cards for 

each cylinder) 



1-10 



110 



Station number 



11-20 FlOoO 



21-30 

31-40 



Note : 



1-10 

11-20 



Mean temperature 

Temperature at inside wall 
minus temperature at outside 
wallo 

Total axial load imposed 
upon cylinder due to all 
causes except pressure „ 
Positive if causing tensiom 
Note' that this is not load 
per unit circumference , but 
total load 0 

There must be Cards 4 for the first and last stations <> 

If no card is submitted for a series of stations, a 

linear interpolation of the temperatures and load is 

made between stations that have cards 0 The cards 

must be in sequence „ Same minimum and maximum number 

of cards as with Card 3 » 

Card Type 5 (No cards or two cards) 
(Used only when Option 2 is used for pressure stresses ») 

F10 o 0 Radius of meridian at junction <, (Ri in Figure 

A2#) Infinity substituted if blank „ 

FlOoO Radius of circumference at junction divided by 

the cosine of the angle between the tangent to 



at this station 
FlOoO Radial tempera- 
ture difference 

FlOoO Axial load 
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the meridian and the axis» (R 2 in Figure A2,) 
Radius of junction substituted if blank « 



Note: There must be Cards 5 after the last Card 4 for 

Cylinders 2 and 3 whenever the pressure Option 2 is 
usedo 

— ___ Skirt 




Figure A2„ Shell, Head, and Skirt » Showing 
R^ and R^o 

Following is an example of the sequence of data cards 
for the input of one problem to the program „ 
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Pressure stress 


Pressure stress 




calculated by 


calculated by 




Option 1 


Option 2 




Card 1 


Card 1 




2 


2 




3 


3 




3 

4 


3 

4 


Cylinder 1 


4 


4 




2 


2 




3 


3 




3 


3 




4 


4 


Cylinder 2 


4 


4 




no card 


5 




2 


2 




3 


3 




3 


3 




4 


4 


Cylinder 3 



4 


4 


no card 


5 



-Denotes the series of Cards 3 and 4 necessary to describe 
the problem „ 

If the problem involves fewer than three cylinders, the 
effects of the cylinders not to be considered may be elimina- 
ted by setting their length equal to a very small number,. 

The program will still give results for these shortened 
cylinders but they might be inaccurate due to round off 
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error and should be ignored . 

The minimum fraction (columns 61-70), Card 1) and 
exponent for forcing convergence (columns 16-20, Card 2) 
are values used in an averaging process to assist in 
obtaining convergence of the numerical solution of the 
fourth order differential equation of a cylinder,, These 
values are named C m; j_ n and m respectively in the main body 
of the thesis o 

The convergence criterion (columns 41-50, Card 2) is 
multiplied by the maximum radial deflection of the cylinder 
and this value is compared with the change in radial de- 
flection at each station for the last two successive 
iterations o If this change for every station is less than 
the maximum deflection times the convergence criterion, 
convergence is deemed to have been attained,, 

The program is stopped by placing a blank card in 
place of Card 1. 
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APPENDIX B 



Listing of Program Skirt and Logical Flow Chart 



List 


of Variables in Main 


Body of Program 


Name in 
Program 


Counterpart 
in Thesis 


Description 


A 


Numerous 


Inputs to continuity and 
statics equations, (Equa- 
tions ( 46 ) in thesis, 
reduced to 6X6 array,) 


AA 


- 




Instruction code for inte- 
gration subroutine "Zigzag” 


AAA 






0 if no divergence exper- 
ienced in problem; 1 if 
divergence experienced , 


AINC 


- 




Number of increments , 


AL 






Coefficient of linear 
thermal expansion , 


AMBAR 


M* 


- 


Junction end moment plus 
the moment necessary to 
compensate for that imposed 
by the radial thermal 
gradient. 


AMO 


Mi 




Moment at junction end. 


AMOO 


M oo 




Externally applied moment 
to junction. 


AV 


A 




Deflection at junction end. 


AX 


V ol> 

M ol> 


V o2> V 

M o2- M o3 


, The end shear and moment 
for each cylinder. 


B 


Numerous 


Inputs to continuity and 
statics equations, (Equa- 
tions ( 46 ) in thesis, 
reduced to 6X6 array,) 



B1 



BM 


M, M* 


Moment in cylinder wall per 
unit circumference. 


BM1 


m l 


Moment in cylinder wall 
caused by thermal, pressure 
and axial loadings. 


BM2 




Moment in cylinder wall 
caused by applying a unit 
shear to the junction end. 


BM3 


M 

rn 


Moment in cylinder wall 
caused by applying a unit 
moment to the junction end. 


BV 


B 


Slope of a straight line 
drawn between radial defies 
tions plotted for junction 
and remote ends. 


CURVE 


X, 


Axial curvature of cylinder 
wall. 


CX 


(k) (x) 


cx = Cc) (x) 


CXI 


- 


CXdX 


CX2 


- 


^ X CXldX 


C 


k 


Pressure necessary to in- 
crease the radius of the 
cylinder one inch. 


Cl 


- 


f 0 X CdX 


C2 


- 


/o X ddX 


E 


E 


Young’s Modulus. 


El 


D 


Axial stiffness of a unit 
width segment of the 
cylinder wall. 


ERROR 


- 


Convergence criterion 
fraction. 


FM 


C 0 
mm 


Minimum fraction in conver- 
gence assistance routine. 
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FRAC 


C 


Fraction used in conver- 
gence assistance routine * 


FX 


f(x) 


Deviation of radial deflec- 
tion from a linear function* 


FXC 


(k)f (x) 


FXC = (FX) (C) 


FXC1 


- 


^ X FXCdX 


FXC 2 




/ Q X FXCldX 


FIX 


f i< x > 


f 0 X GXdX 


F2X 


f 2 (x) 


Normalized FIX* 


GX 


g(x) 


f Q x (BM/D)dx 


H 


t 


Thickness of cylinder wall* 


ICYL 


- 


Cylinder number* 


IEXP 


m 


Exponent in convergence 
assistance routine* 


IIF 


- 


Printout option* 


INC 


n-1 


Number of equal increments 
into which the cylinder is 
divided* 


INC2 


- 


(2) (INC) + 1 


IND 


- 


Interval of printout of 
input data* 


IPROB 


- 


Problem number* 


IWILL 


Option 1 or 2 


Pressure stress calculation 
opti on* 


JJ 


- 


Code for type of calculation 
presently in progress* 


K1 


— 


First station in a sequence 
in which the loadings were 
not read into the computer* 
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K2 




First station in a sequence 
in which the loadings were 
read into the computer,, 


K3 




(K2 - K1 + 2)/2 


L 


- 


Station number where maxi- 
mum F2X occurs. 


LL 


- 


Station number where maxi- 
mum deflection occurs „ 


NAA 


- 


Number of Iterations „ 


PSI 


P 


Internal pressure * 


PSI1 


P 


Internal pressure * 


R 


a 


Radius of cylinder,. 


R1 


r i 


Meridional radius of curva- 
ture of at junction* 


R2 


R 2 


Radius of curvature of line 
perpendicular to meridian 
at junction* 


SIG1 




Longitudinal membrane stress 
due to pressure* 


SIG2 




Circumferential membrane 
stress due to pressure* 


SL 


ds 

dx 


Axial slope of cylinder wall 


SLOT 


- 


Value of X at end of slots* 


SLOT1 


- 


Value of X at beginning of 

SlOtS a 


ST 


<*E4T 

2 


Stress at surface of wall of 
an infinite cylinder sub- 
jected to a radial tempera- 
ture gradient* 


STCI 


(T . 
z inner 


Circumferential stress on 
inner wall surface* 
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STCO 


cr 

z outer 


Circumferential stress on 
on outer wall surface „ 


STLI 


cr • 

^x inner 


Longitudinal stress on 
inner wall surface „ 


STLO 


outer 


Longitudinal stress on 
outer wall surface., 


TL 


T 


Mean temperature „ 


TLINC 


- 


Increment cf linear inter- 
polation for TL, 


TP 


— 


Number of 53 word blocks 
needed to store one para- 
meter of cylinder,. 


TPI 




Total number of 53 word 
blocks necessary to store 
on tape the intermediate 
values of one problem „ 


TR 


A T 


Difference between inner 
and outer wall temperature „ 


TRINC 


- 


Increment of linear inter- 
polation for TRo 


UM 


M 


Poisson’s ratio,, 


V 


V 


Radially outward directed 
load per unit circumference 


VBAR 


- 


Constant shear imposed at 
junction end„ 


VBAR1 


- 


^ X VBARdX 


VO 


v i 


Shear at junction end., 


VOO 


V 

oo 


Externally applied shear 
force to junction c 


VI 


t l 


Shear in cylinder wall 
caused by thermal, axial, 
and pressure loadings » 
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V2 


V v 


Shear in cylinder wall 
caused by applying a unit 
shear to the junction end. 


V3 


V m 


Shear in cylinder wall 
caused by applying a unit 
moment to the junction end. 


WT 


(P*)( 27 ra) 


Total axial load caused by 
all influences other than 
pressure , 


WTINC 


- 


Increment of linear inter- 
polation for WT, 


W 


Q+££(A+Bx+f(x) ) 
a- 


Distributed axial loads 
applied to cylinder includ- 
ing simulated thermal loads 


W1 


- 


f X WdX 


W2 


- 


^ X V/ldX 


X 


X 


Axial dimension. 


XDEL 


- 


Length of interval between 
stations. 


XMAX 


x n 


Value of X at remote end. 


XMIN 


X 1 


Value of X at junction. 


Y 


s 


Radial deflection. 


Y1 


S L 


Radial deflection of the 
cylinder wall caused by 
thermal, pressure, and axia 
loadings. 


Y2 


s v 


Radial deflection of the 
cylinder wall caused by 
applying a unit shear to 
the junction end. 
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Y3 s m 

m 


Radial deflection of the 
cylinder wall caused by 
applying a unit moment to 
the Junction end* 



The following variables are those used in Subroutine Zigzag* 



Variable 


Description 


A 


Instruction code for integrationo 


B 


Interval length * 


D 


Integrand* 


E 


Integral* 


M 


(2) (number of increments) + 1. 



The following variables are those used in Subroutine Answer * 



A 


Equivalent stress on outer 
wall surface* 


B 


Equivalent stress on inner 
wall surface* 


BM 


Moment » 


C 


Outer longitudinal stress* 


D 


Inner longitudinal stress* 


E 


Outer circumferential stress* 


F 


Maximum equivalent stress* 


G 


Inner circumferential stress* 


KK 


Denotes whether maximum equivalent 
stress is on inner or outer surface 
of wall* 


KKK 


Station number of location of 
maximum equivalent stress* 


L 


Cylinder number* 
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M 



Number of increments . 



N 


Printout code. 


NN 


Problem number. 


V 


Shear. 


X 


Longitudinal dimension. 



The symbols used in the following logical flow chart are 
from the '’Flowcharting Template,” Form X20-S020, of Internal 
tional Business Machines Corporation, and are described below. 
The numbers in the locator (circle) and connector (inverted 
house) correspond to the sequence numbers at the extreme 
right of the program listing. 



Symbol 


Use 


Oval Terminal 


The beginning of end of a program. 


Circle Locator 


The location of the first statement 
of the following operation in the 
program. 


Trapezoid Input/Output 


Any function of an input/output devic 


Inverted Connector 

House 


An entry from, or an exit to, another 
part of the program flowchart. 


Diamond Decision 


A decision is made for routing to 
different branches of the program. 


Rectangle Processing 


A group of instruction performing 
a processing function. 



Subroutine GAUSS2 is not charted or listed because it is 
a standard library program for solution of simultaneous 
linear equations. The subroutine was written by C. B. Bailey 
and is listed in ”F2 UTEX LINEQN,” CO OP Manual for CDC 1604, 
December, I960. 
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SKIRT (continued) 



= 1 






! = 1 

j 



CYLINDER 
NUMBER 



r 3 




ENTER INPUTS 
TO CONT. AND 
STAT. EQNS. 
WITH LOADS 
APPLIED TO 
CYLINDER 1 



SAME INPUTS BUT 
LOADS APPLIED 
TO CYLINDER 2 



AME INPUTS BUT 
iOADS APPLIED 
[TO CYLINDER 3 



SAME INPUTS RUl] 
UNIT SHEAR 
APPLIED TO 
CYLINDER 1 



-= 2 _ 



SAME INPUTS 
UNIT SHEAR 
A P -LIED TO 
CYLINDER 2 



BUT 



SAME INPUTS 
UNIT SHEAR 
APPLIED TO 
CYLINDER 3 



BUT 



CORRECT Y, 
AND BM. 
CALCULATE ST. 
w r o. 



WRITE TAPE 
\ 1 ST,E. 



SAME INPUTS BUI 
UNIT MOMENT 
APPLIED TO 
CYLINDER 1 




SAME INPUTS BUI 
UNIT MOMENT 
APPLIED TO 
CYLINDER 2 



SAME INPUTS BUT 
UNIT MOMENT 
APPLIED TO 
CYLINDER 3 
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SUBROUTINE 

ZIGZAG 



SUBROUTINE 

READ 



SUBROUTINE 

ANSWER 





i. 
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APPENDIX C 



Sample Problems 

Elementary problems were set up to show the capabili- 
ties of the program. No pretense is made that the values 
used are representative loadings. The problems are des- 
cribed below. No significance should be attached to the 
numbering system of the problems. A summary of the equi- 
valent stress in each cylinder of each problem is given in 
Table 1 and sample printouts of problems 20 ,21 ,22, SO, and 
90 along with the input cards for these problems, are given 
after page C6. 

Problem 10 - Three identical cylinders, 16 inches long, 50 
inch radius, 0.5 inch wall thickness. Young *s Modulus 
3X10'' psi, Poisson's Ratio 0.3, thermal coefficient of 
expansion 7.7X1Q~6 inches/inch ®F, number of increments 
25 o 100 psig internal pressure in Cylinders 2 and 3, no 
thermal or mechnical loadings other than pressure. The 
method of solution used was Option 1 (See Appendix A). 
Printout option 3. 

Problem 20 - Same as problem 10 except Option 2 used for 
calculation. In cylinder 2, R^ - 50 inches, R 2 “ 50 
inches; i.e., hemispherical head. In Cylinder 3, Rp is 
infinite, R 2 ~ 50 inches. 

Problem 21 - Same as problem 20 but printout option 2 is 



Cl 



used* 

Problem 22 - Same as problem 21 but printout option 1 
Problem 30 - Same as problem 20 with following thermal 
loading added" 

Cylinder 1 Medial wall temperature = T - 300 - 
40 Q(l~e^ 025x ) °F U 



Cylinders 2 and 3 T = 300 + 200(l~e~ o025x ) ®F» 

Problem 40 = Same as problem 30 with the following thermal 
loading added; 

Cylinder 1 Axial temperature difference -^T - 0 for all x. 
Cylinders 2 and 3 4T = ( 25 ) C T i,/Ti_) 

where T. is the median temperature at station i and T% is 

the median temperature at the junction. 

Problem 50 - Same as problem 40 with an axial load added; 

Cylinder 1 Axial load » P* = =5000 - 47 ^5x Lb. 

Cylinder 2 P* = 3 500=47 . 5x Lb. 

Cylinder 3 P* = =1500+47 »5x Lb. 

Problem 60 - Same as problem 50 but the thickness of 

Cylinder,^ is progressively reduced. 

Thickness of Cylinder 3 
0.0f x^2.56 .5 Inches 



2.56£x*5o76 

5.76£x*3c96 
8.96^x^12.16 
1 2 . 16-xf 16 . 0 



.45 Inches 
.40 Inches 
.35 Inches 
.30 Inches 



C2 



Problem 70 - Same as problem 60 but Cylinder 1 (skirt) is 

•slotted from x = 0 to x = 10 inches* 

Problem SO - Same as problem 70 but externally applied 

moment and shear are addedo 

External moment = M ~ 500 lUrlb 

oo in 



External shear ~ V — 500 

oo m 

Problem 90 = Same as problem 20 but length of cylinder 1 is 



reduced to 10“^ inches, 



effectively eliminating it from 



the problem o 
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Equivalent Stresses in Sample Problems 
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Discussion of Results 

Since the sample problems are somewhat complex and are 
used mainly to display the capabilities of the program, the 
purpose of the following discussion is mainly to point out 
those results which the writer felt most interesting and 
to offer what seem to be plausible explanations of these 
re suit So 

Problems 20 and 90 have no loadings but pressure and 
differ only in that 20 has a skirt and 90 does not In 
comparing the stresses in Table 1 of Cylinders 2 and 3 > it 
will be noted that the addition of the skirt did not 
grossly affect the stresses or their locationso The 
stresses did increase somewhat due to the moment exerted 
by the skirt upon the junction, but this increase was 
partially balanced by the skirt*s constraint on the outward 
movement of the junction * This may be seen by comparing 
the circumferential stresses of Cylinders 2 and 3 of the 
two problems o The results for Cylinder 1 of problem 90, 
should be neglected » 

Note that the application of an external moment and 
shear on problem SO considerably changes the equivalent 
stress at the remote end of Cylinder 3 from what it was in 
problem ?0o. Since the length of the cylinder was supposedly 
chosen such that a force applied at the junction should not 
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affect the remote end, this change in stress might be 
puzzling* However the junction end moment in Cylinder 3 
of problem 70 was “8$ in-lb* and in problem 80 was 559 
in-lb o, an increase of over 500%* Also the junction end 
shear increases from 59 Ib/in to 471 lb/in, an increase of 
about 700%* Since, due to the sign conventions for Cylinder 
3, (See Figure 5 of text) a moment and shear of opposite 
signs have a cumulative effect on the moment, this 14 
percent change in stress at the end is not unreasonable* 

From the printout of problem 80 it will be noted that 
at station 25 of Cylinder 3 the stresses seem to change 
radically from those of the previous station* The probable 
reason for this is that in this particular cylinder a large 
nonlinearity is imposed upon the whole length of the 
cylinder by the decrease in wall thickness. Such a non- 
linearity will create internal moments and shears in the 
wall* However we are still forcing the shear and curvature 
to be zero at the end (Station 26) and therefore a large 
compensation occurs near the end* For accurate results 
near this point, a longer cylinder must be used with more 
intervals* 
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DATA CAROS 
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